Spin Dynamics of La2Cu04 and the Two-Dimensional Heisenberg Model 

Anders W. Sandvik(°)'(''' and Douglas J. Scalapino 
Department of Physics, University of California, Santa Barbara, CA 93106 

(October 28, 1994) 



> 

o 



(N 

O 

o 



o 
o 
I 

Oh 



The spin-lattice relaxation rate 1/Ti and the spin echo 
decay rate I/T2G for the 2D Heisenberg model are calculated 
using quantum Monte Carlo and maximum entropy analytic 
continuation. The results are compared to recent experiments 
on La2Cu04, as well as predictions based on the non-linear 
o"-model. 

The Cu02 planes of the undoped high-Tc cuprates are 
good physical realizations of the two-dimensional (2D) 
antiferromagnetic Heisenberg model. [1] The mapping of 
this lattice model onto the non-linear cr-model (nlcrm) 
in 2-1- 1-D has led to detailed predictions for various ex- 
perimentally measurable quantities. [2-4] For T < QOOK 
(above the 3D ordering temperature) the correlation 
length of La2Cu04 grows exponentially as the temper- 
ature is lowered. [5] The behavior is in close agreement 
with Quantum Monte Carlo results for the 2D Heisenberg 
model with a nearest-neighbor coupling J « 15007i', [6] 
and corresponds to the the nlcrm in the low-temperature 
"renormalized classical" (RC) regime. [2] It was recently 
suggested [3,4,7,8] that the high-temperature behavior of 
the cuprates corresponds to the "quantum critical" (QC) 
regime of the nlcrm, where the leading temperature de- 
pendence of the inverse correlation length is linear. Ex- 
perimental evidence supporting this scenario has been 
provided by Imai et ai, who measured the spin-lattice 
relaxation rate 1/Ti [9] and the gaussian component of 
the spin-echo decay rate I/T2G [10] at temperatures as 
high as T = 900K. In particular, Imai et al. found that 
1/Ti and the ratio T1T/T2G were both temperature in- 
dependent at high temperatures, as predicted for the QC 
regime. [8] These experiments were recently repeated by 
Matsumura et al.. [11] Their results for 1/Ti are almost 
identical to the earlier ones, but for I/T^g the temper- 
ature dependence obtained is different at high temper- 
ature, causing TT1/T2G to be temperature dependent, 
in disagreement with the QC scenario. In addition to 
this discrepancy, an open question is the reason for the 
absence of the minimum in 1/Ti at T « 750K, theoreti- 
cally predicted by Chakravarty and Orbach. [12] In order 
to settle these questions we have calculated both 1/Ti 
and I/T2G for the 2D Heisenberg model using quantum 
Monte Carlo simulation and the maximum entropy ana- 
lytic continuation method. This enables us to compare 
directly the spin dynamics of the Heisenberg model and 
La2Cu04, as well as to assess rigorously the accuracy 
of the predictions based on the nlcrm. Lower tempera- 
tures can be reached than with high-temperature series 
expansions [15] and calculations on small clusters, [16] 



and the approximations necessary with these methods 
can be avoided. 

We find that I/T2G for the 2D Heisenberg model de- 
creases faster than the rate measured for La2Cu04 by 
Imai et al above 750K. The temperature dependence is 
~ T"2 for 0.45 < T/J < 1, in disagreement with the 
QC prediction ~ T^^. Taking into account the temper- 
ature dependence of the spin wave velocity [17] slightly 
improves on the agreement. For 1/Ti our results are in 
good agreement with the experiments. We also find that 
while 1/Ti exhibits a minimum for a local contact hy- 
perfine coupling of the type used by Chakravarty and 
Orbach [12], this minimum is absent when the experi- 
mentally known on-site and near-neighbor interaction is 
used. 

We study the standard 2D Heisenberg hamiltonian 



H 



J 2_^ /__/ ^i ' ^t+6 



(1) 



where Si is a spin-i operator, and S runs over the near- 
est neighbors of site i. The spin-lattice relaxation rate 
and the spin-echo decay rate for a given nucleus provide 
information on the spin susceptibility through the direct 
and transferred hyperfine couplings of the nuclear spin to 
surrounding electronic spins. For a ^^Cu nuclear spin Jo 
at site 0, the coupling to the electronic spins Si is given 
by the hyperfine hamiltonian [18,19] 

''H = A^mSS + I'oS'o) + A\\ISS^ +bY,Io-Ss. (2) 



The constants ^x, ^|| and B are known from Knight shift 
measurements. 

With the external field in the direction a, the NMR 
spin-lattice relaxation rate is given by [20] 
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where a' denotes the two axes perpendicular to a, and 
A^ is the fourier transform of the a' component of the 
hyperfine coupling. The dynamic structure factor S{q,uj) 
is related to the imaginary part of the spin susceptibil- 
ity; S{q,uj) = x"(9,i^)/(l — e~^"). Since the resonance 
frequency coff is small compared to J, 1/Ti effectively 
measures S{q,uj —>■ 0), averaged with the hyperfine form 
factor |A™ 1^. In terms of the inverse fourier transform 
Smn = S{mx + ny,uj —^ 0) of S(q,uj), 1/Ti with the ex- 
ternal field perpenducular to the Cu02 planes is given 
by 



+ 16A±BSw + WB'^Sn + 8B'^S2o. 



(4) 



This is the rate measured in the experiments by Imai et 
al. and Matsumura et al.. Smn{'^) can be obtained from 
the imaginary-time correlation function 

Cmr.{T) = (5^,+„^(r)5onO)), (5) 

where Sp{T) = e'^^ Spe~'^^ , by inverting the relation 



(-'mn\T) 
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(6) 



Here Cmn{T) is computed using quantum Monte Carlo 
simulation, and the inversion of (6) is carried out with 
the maximum entropy method. [13,14] 

The rate I/T2G is related to the interactions between 
the nuclear spins. The coupling (2) leads to an indi- 
rect nuclear spin-spin interaction, which dominates the 
direct dipole-dipole interactions. For the external mag- 
netic field applied perpendicular to the CUO2 planes, 
Pennington and SHchter derived the following expression 
for 1/T2G- [21] 
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(7) 



Here Jz{xi) is the z-component of the induced interaction 
at distance Xi, given by 

Mxi) = AiiF,{xi) + Bj2FM+s), (8) 



with 



where x(xi) is the static response at separation f= Xt, 
given by the Kubo formula 



Fz{xi) = -^ (A||x(f,) +bY, x{S!r+s)) , (9) 



X{x 



) = jdT{st{T)s^,m. 



(10) 



The factor 0.69 in (7) is the natural abundance of the 
^^Cu isotope. 

We have used a recently improved version of the Hand- 
scomb quantum Monte Carlo technique [22] to calcu- 
late the necessary correlation functions. Unlike standard 
methods, [23] this technique is free from the systematical 
errors associated with the Trotter break-up. For the an- 
alytic continuation of the imaginary-time data necessary 
to obtain 1/Ti, we have implemented the so called "clas- 
sic" maximum entropy procedure as described in a recent 



paper by Jarrell and Gubernatis. [14] We have studied 
systems of iV = 64 x 64 spins with periodic boundary 
conditions, at temperatures T / J = 0.25 — 1.0. At these 
temperatures the correlation length is smaller than the 
lattice size, and there are virtually no finite size eff'ects. 
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FIG. 1. Monte Carlo results for I/T2G (solid cir- 
cles) and experimental results by Imai et al. [10] (open 
squares) and Matsumura et al. [11] (open circles). The 
solid line in the main figure is of the form ~ T^^. In- 
set: Same as the main figure with the QC prediction by 
Chubukov et al. [24] (solid curve) and including the tem- 
perature dependence of the spin wave velocity (dashed 
curve). 

The calculation of I/T2G is straight-forward, as it in- 
volves only the static susceptibility (10). We use the 



relation A 



-4:B, experimentally known to hold quite 



accurately. [18,19] We are then left with J and B as fit- 
ting parameters, that can be checked against other exper- 
iments. The best agreement with the experimental data 
for I/T2G is obtained with J = 1580K and B = 3.4x10"'' 
eV (« 37 kOe/fis), both consistent with other estimates. 
[5,6,19] The Monte Carlo results with these parameters 
are shown in Figure 1, along with the experimental data. 
Although the over-all ageement is good, a notable fea- 
ture is that for T > 750K the data of Imai et al. [10] 
is flatter than both the Monte Carlo results and those of 
Matsumura et al.. [11] This flatness cannot be reproduced 
for the Heisenberg model with any reasonable values of 
J and B and, if correct, must be associated with physics 
not described by this model alone. On the other hand, 
the data of Matsumura et al. is well reproduced at high 
temperatures, but exhibits slightly more curvature than 
the Monte Carlo results and the data of Imai et al. at 
intermediate temperatures. Fig. 1 also shows the theo- 
retical form derived by Chubukov et al. [24] for the QC 
regime, which for the hyperflne couplings used here be- 
comes 



1 

T2G 



0.491 ^(T) X 10" s" 



(11) 



where the QC correlation length is given by [3,4] 



$, = c/(1.038 T) (QC regime), (12) high-temperature limit. The RC expression is [12,4] 



and the spinwave velocity c « 1.68. [27] As noted by 
Chubukov et al., [24] the over-all magnitude of I/T^g is 
well reproduced with this formula, but the slope is not. 
Actually, the Monte Carlo results for I/T2G in the regime 
0.45 < T < 1 is very well described by a 1/T^ behavior; 
a quite significant deviation from (11). Elstner et al. [17] 
recently pointed out that the leading lattice corrections 
to the nlcrm can be taken into account via a tempera- 
ture dependent spinwave velocity. The velocity calcu- 
lated from Monte Carlo results for the static structure 
factor and the static susceptibility agrees well [25] with 
the high-temperature series expansion results by Elstner 
et al., [17] and when used in Eq. (12) slightly improves 
the agreement with the Monte Carlo results for I/T2G. 

We now turn to the calculation of 1/Ti, which is more 
complicated as it relies on a numerical analytic contin- 
uation of imaginary time correlation functions. For the 
local correlation function Coo(t) the relative statistical 
errors in our data is typically as low as 10"'', and the 
continuation of this quantity is relatively stable. For Cio, 
Cii and C20 the relative errors are typically on the or- 
der of lO^'^ and an accurate determination of 1/Ti using 
the full extended hyperfine couplings is therefore more 
difficult than with a strictly local one. 



0.2 - 



0.15 



0.1 - 



0.05 




FIG. 2. Maximum entropy results for Sqo vs T (sofid 
circles) compared to the RC form (13) and the high- 
temperature form (15) (sofid curves). The dashed curve 
is the QC form (16). 

In Figure 2, the to —>■ limit of ^oo is graphed ver- 
sus the temperature. For a strictly local coupling, this 
quantity is proportional to 1/Ti. Repeating the ana- 
lytic continuation procedure for different subsets of the 
Monte Carlo data, we estimate the statistical errors to 
be approximately 5-10% (any bias due to the maximum 
entropy procedure itself is of course not captured this 
way). A broad minimum around T/J = 0.5 is observed, 
in good agreement with the prediction by Chakravarty 
and Orbach, [12] who deduced this feature by contrast- 
ing the behavior of Soo{uj —>■ 0) in the RC regime and the 
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where the correlation length is given by [2,26] 
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2tvp, 
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T 

4tvp, 



J^^Ps/T 



(RC regime). (14) 



The spin-stiffness ps « 0.18 [27,28] and the ordered mo- 
ment is No « 0.31. [28] The constant A has not been cal- 
culated rigorously, but an estimate based on fitting the 
nlcrm scaling forms to numerical results is XNq = 0.61. 
[12,4] The rather poor agreement with our result for 5*00 
shown in Figure 1 indicates that this value is too large. It 
should be noted, however, that even the lowest temper- 
atures studied here correspond to the cross-over regime 
to RC behavior, [17] and perfect agreement with the RC 
expression cannot be expected. The high-temperature 
form is [12] 



5o^o^(c. ^ 0) 



(0F/8)e 



-{l/2Tf(l + l/4T) 



7l + l/4T + 0[T-2]) 



(15) 



and deviates from the Monte Carlo data by 25% at T/J = 
1. Clearly terms of order T^^ in Eq. (15) might be 
important for T < 1. In the regime 0.4 < T/J < 0.7 
S'oo(w -^ 0) appears to be rather flat, as predicted for 
the QC regime. The expression derived by Chubukov et 
al. is, [4] 



Sg'^iLO ^ 0) 



M(JL 

Ps 27rps 
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(16) 



where the the 3D classical Heisenberg exponent rj « 0.03. 
Deep inside the QC regime Ri is a constant, for which 
Chubukov et al estimated Ri « 0.22 (there are certain 
complications in estimating i?i). [4] Eq. (16) with this 
value of Ri describes the behavior in the intermediate 
temperature regime quite well. 




-100 



600 



800 



1000 



1200 



FIG. 3. Maximum entropy results for 1/Ti vs T (solid 
circles) compared to the experimental results by Imai et 
al. [9] (open squares) and Matsumura et al. [11] (open 
circles). 



The minimum in 1/Ti has not been observed experi- 
mentally. [9,11] In Figure 3 we show results obtained with 
the full hyperfine coupling (2), using the same values of 
J and B as in the fit to I/T2G in Figure 1. For A±/B we 
take the experimental value 0.84. [18,19] The statistical 
errors are rather large, as discussed above, but a clear dif- 
ference from the temperature dependence of Figure 2 can 
be noted, and the agreement with the results by Imai et 
al. [9] and Matsumura et al. [11] is reasonably good ( an 
even better agreement is obtained with a slightly smaller 
A±/B). In particular, 1/Ti is temperature independent 
at high temperatures and the minimum found above for 
Soo{uj —>■ 0) a.t T » 800K is not reproduced. 

To summarize our results, we note that I/T2G for the 
Heisenberg model agrees well with the experimental re- 
sults for La2Cu04 by Imai et al. [9,10] and Matsumura 
et al. [11] for T < 750K. However, for T > 750K 
the Heisenberg result decays considerably faster than the 
data by Imai et al.. The temperature dependence is 
close to 1/T^ in a wide temperature regime. For 1/Ti 
our results are also in reasonable agreement with the 
experiments and, in particular, is almost temperature- 
independent at high temperatures, as predicted for the 
QC regime. [4] As a consequence of the bahavior of 
I/T2G, the ratio T1T/T2G is not constant at high tem- 
peratures, in contradiction to the QC prediction. [8] It 
appears that the impHcations of the high-temperature 
results for I/T^g by Imai et al. have to be reconsidered. 

In a recent paper, Elstner et al. [17] found that pure 
QC behavior for the correlation length of the 2D Heisen- 
berg model can be observed only above T/J « 0.6, and 
that lattice effects also become important at these tem- 
peratures. In the temperature regime 0.4 < T/J < 0.6 
there is a cross-over to RC behavior. In view of these re- 
sults the poor agreement with QC scaling for I/T2G found 
here is perhaps not surprising. On the other hand, 1/Ti 
does exhibit QC behavior for T > 700K. It should also 
be noted that the uniform susceptibiHty exhibits strik- 
ingly accurate QC behavior at lower temperatures. [3,4] 
Hence, it is clear that effects of the proximity to the criti- 
cal point are manifest in the 2D Heisenberg model. How- 
ever, the size and location of the regime where QC be- 
havior can be observed depends strongly on the quantity 
considered, and there does not seem to exist a tempera- 
ture regime where the QC scaHng formulas can be appHed 
universally. 
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